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Abstract 

We study the evolution of branching trees embedded in Euchdean spaces with suppressed branch- 
ing of spatially close nodes. This cooperative branching process accounts for the effect of over- 
crowding of nodes in the embedding space and mimics the evolution of life processes (the so-called 
"tree of life" ) in which a new level of complexity emerges as a short transition followed by a long 
period of gradual evolution or even complete extinction. We consider the models of branching trees 
in which each new node can produce up to two twigs within a unit distance from the node in the 
Euclidean space, but this branching is suppressed if the newborn node is closer than at distance 
a from one of the previous generation nodes. This results in an explosive (exponential) growth in 
the initial period, and, after some crossover time tx ~ ln(l/a) for small a, in a slow (power-law) 
growth. This special point is also a transition from "small" to "large words" in terms of network 
science. We show that if the space is restricted, then this evolution may end by extinction. 
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I. INTRODUCTION 



"The evolution of life is, obviously, a nonuniform process" [l| . For biological evolution, 
this means that new types of biological objects emerge abruptly with subsequent gradual 
evolution. This evolutionary process can be schematically depicted as a tree ("the tree of 
,|2|), where branches, are, for example, different species. Importantly, the growth of 



life 



this tree is complicated by interaction and competition between species. In this work we 
discuss one of the simplest models of growing trees which can mimic this process. Growing 
trees naturally represent a wide range of real- life processes and phenomena js 10|. 

n 

The Galton- Watson branching process [111] provides a simple example of a growing tree 
with non-interacting nodes and so uncorrelated branching. Interacting branching processes 



are much more interesting and difficult for analysis [12|]. In this work we study evolving trees 
which evolution is influenced by interaction between some of existing nodes, for example, 
nodes of the previous generation. We assume that the growing tree is embedded in some 
metric space and assume that spatially close nodes of the previous generation suppress 
mutually their ability to born new nodes. In other words, overcrowding of nodes in the 
embedding space suppresses their "fertility". This kind of interaction (competition), leading 
to suppression of branching, emerges if there is no sufficient space, no niches for the new 
nodes and branches (species). 

Due to the embedding space, we can introduce distance between two nodes other then 
the shortest path internode distance for this tree. For the sake of simplicity, we consider 
a D-dimensional Euclidean space, although the results do not depend qualitatively on D. 
Networks embedded in metric spaces and their evolution already attracted much attention 
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18| . In this paper we are particularly interested in a transition (actually, crossover) 
between different regimes of the network growth, namely, explosive (exponential) evolution 
and gradual (power-law) one. Here the evolution of the network is characterized by the 
variation of the number of its nodes (which corresponds to biological diversity, for example). 
We find the position of this transition and express it in terms of a single model parameter. 
This transition coincides with crossover from a "small word" to "large word" network ar- 



chitectures 



19| , where small words show a logarithmic dependence of network diameters on 



their sizes (total numbers of nodes) and large words show a power-law dependence 20|, |21 |. 



One should emphasize a principal difference from the previous studies of this crossover. 




FIG. 1. (Color online) The scheme of the network growth on a plane. The black node shows the 
root, the nodes of the first and the second generations are dark and light grey, respectively. The 
furthermost left node attempts to born two children. The first attempt (a) is abandoned because 
of the nearby second generation node. The second attempt (b) is successful since the new node 
has no second generation node within radius a from it. This results in the network (c). 

In Ref. [l^, the crossover was controlled by a model parameter, while in the present study 
the small-world and large-world architectures are realized on different stages of the network 
evolution. In addition we find how the spatial distribution of nodes evolves. We also 
consider the evolving trees embedded in restricted areas of metric spaces, and investigate 
the possibility of complete extinction. The term extinction for these trees implies the end 
of evolution and the absence of new generations. 

The paper is organized as follows. In Sec. II we present the model of an evolving tree 
with interacting nodes and obtain simple estimates for its growth. In Sec. Ill we describe 
the network size evolution for two particular models of these trees. In Sec. IV we consider 
the trees embedded in a bounded space and extinction. In Sec. V we describe the evolution 
of the node spatial distribution. Finally, in Sec. VI, we summarize our results. 



II. THE MODEL 

The model, which we use, is formulated as follows (see Fig. [1] showing schematically the 
grows of the tree embedded in a two-dimensional space). The growth of the tree starts from 
a root node. At each time step, each of the nodes of the tree attempts to emit two leaves 
(leaf is a link with a new node), so at each time step a new generation of nodes is born. The 
network is embedded in a D-dimensional Euclidean space, and the root has zero coordinates. 
At each time step, make the following: 
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(1) Choose uniformly at random a node i (coordinates Xj) from the previous generation 
and make an attempt to create its leaf with a new node at the point Xj + Aj. Here the 
random vector Aj is uniformly distributed within — l<A^j<l, — l<Aj^j<l, .... If 
among the nodes of the previous generation (excluding the parent node i) and among 
the nodes already created at this time step, no nodes are closer than at distance a from 
the point, Xj + Aj, then create the leaf. If such nodes exist, abandon this attempt. 
Make the next attempt to create the second leaf from this node using the same rules. 

(2) From the rest nodes of the previous generation, choose uniformly at random nodes one 
by one and repeat (1) again and again until all the nodes of the previous generation 
will be updated. 

We will also consider a variation of this model, in which for each attempted node birth, 
closeness to all existing nodes should be checked and not only to the previous generation 
nodes. Importantly, in both these models, existing nodes and links never disappear. 

Figure [2] shows the result of simulation of this model for D = 1 and sufficiently small a, 
namely, the evolution of the number of nodes N{t) of generation t which plays the role of 
time. Initially, grows exponentially, N = 2^. One can see that after certain crossover time 
tj., the network growth is slower than exponential. For an arbitrary dimension D, one can 
easily estimate = const (t/a)^ at large t. To obtain this estimate, we assume that nodes 
of generation t are within a hypersphere which radius grows with a constant rate of the 
order of 1 (the rate is actually smaller than 1). This average rate of expansion is explained 
by the fact that children in this tree are born within unit distance from their parent nodes. 
Since the neighboring nodes cannot be closer than at distance a, we obtain ~ t^/a^. 

Note that if the parameter a is sufficiently large, does not grow at all. If a is, say 2, 
A^ = 1 for any t, and our tree is a chain. 

From A^(t < t^) = 2* and N{t ^ t^) ~ (t/a)^, we have very roughly 




(1) 



which leads to the estimate 




Time 



FIG. 2. (Color online) The number of new nodes in the tree as a function of time, D = 1 and 
a = 0.1. 

at small a. In Sec. El we will demonstrate that this simple estimate describes well the results 
of our simulations. 

III. NETWORK SIZE EVOLUTION 

Data similar to Fig. |2] are shown on the normal- log plot, Fig. [3]^a), for a few values of 
a {D = 1). The straight line in the figure is the dependence 2*, and the crossover from 
the exponential to a slower growth is clearly seen. Figure EJ^a) was obtained from the 
model formulated in Sec. [Tll in which the previous generation nodes affect the branching 
process. We performed similar simulations for the model, in which newborn nodes cannot 
be closer than at distance a from any of existing nodes (apart of their parents). The results 
of the simulations (the evolution of the number of nodes of generation t) are shown in 
Fig. [31(b). In contrast to Fig. El^a), in the network in which all nodes influence branching, 
the number of nodes of generation t approaches a constant value Nmax{ci) at large t. One 
can easily obtain this plato using an estimate similar to that from the previous section, 
Ntot — const (t/a)^. The only difference is that now Ntot in that estimate is the total 
number of nodes in the network, and so for the number of nodes of generation t, we have 
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FIG. 3. (Color online) Evolution of the number of new nodes in the networks for different values 
of the parameter a {D = 1). The data were obtained after averaging over 100 samples, (a) The 
trees evolve according to the rules introduced in Sec. [Til i-e., only the previous generation nodes 
influence the branching process, (b) The trees in which newborns cannot be closer than at distance 
a from any of existing nodes (apart of their parents), i.e. the branching process is influenced by 
all existing nodes. 

Nt = dNtot{t)/dt ~ Dt^~^ /a^ , and, in particular, Nf ~ 1/a at D = 1. The results of 
simulations for this model, which give Nmax ~ OA/a, see Fig. HJ agree with this simple 
estimate. 

The crossover time between two regimes of the network evolution is obtained in Fig. [5] 
for the growing tree model {D = 1) from the previous section. Note that the result, tx{a) = 
-0.34 + 1.46 ln(l/a), agrees well with Eq. ([2]), U ~ (£>/ ln2) ln(l/a), since l/ln2 = 1.44.... 
Clearly, the time t is of the order of the diameter d of this tree (the maximum separation 
between two nodes in a network). So we have the logarithmic dependence of the diameter d 
on the total number Nfot of nodes in these trees for t <^tx, and the power-law dependence 
d{Ntot) for t tx, which corresponds, respectively, to the small- world and large- world 
network architectures. 

IV. SPATIAL RESTRICTION 

It is natural to introduce a spatial restriction into the models. If our network is embedded 
in an infinite Euclidean space, the evolution is actually determined by the only parameter a 
(recall that we set the scale of displacements of children nodes from their parents to 1, i.e., 
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FIG. 4. (Color online) Log-log plot of N^aax versus a obtained by simulating the model of trees 
in which newborn nodes cannot be closer than at distance a from any of existing nodes (apart of 
their parents). The straight line has slope -1. 



this is the unit distance in this problem). If the area of the space, in which the network is 
embedded, is restricted, then the model has an extra parameter, namely the linear size of 
the area, L. For simplicity, we assume that this area does not change with time. 

Let the Euclidean coordinates Xj of all nodes in the network be within the area — L < 
Xi < L, —L < Hi < L, . . . . In our simulations we use periodical boundary conditions, but, 
in principle, this is not necessary. If L is finite, then one may expect that the size of the tree 
will finally approach some limiting value. The network has even a chance to extinct if at 

aiology, the smaller 



22|. 



some moment all its nodes occur in one small area. In, e.g., population 
a population, the more susceptible it is to extinction by various causes [2 

Figure [6] demonstrates an example of the evolution of the network, which in this case 
has a = 0.1 and L = 1. The network rapidly enters the fluctuation regime, in which 
Nf fluctuates around a mean value Nmax, and extincts before 900 time steps. After that 
we again introduced a root and restarted the process. The picture which we observe agrees 
with traditional views on extinction processes which show "relatively long periods of stability 
alternating with short-lived extinction events" (D. M. Raup) 23(|. This kind of extinction 
may occur in branching annihilating random walks and other related processes studied in 
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FIG. 5. (Color online) Finding the crossover time tx from simulation data for the model from the 
previous section. N = Nt is the number of the t-generation nodes. The dependencies versus 
t — tx for different values of a and D = 1 (a) collapse into a single curve for the crossover times tx{a) 
shown on panel (b). Fitting gives tx{a) = —0.34 + 1.461n(l/a), which is consistent with Eq. ([2]). 



Refs. 2J, |25|. In other models of biological evolution, extinction may require external factors 



or an environmental stress 26'1 or an internal mechanism, such as a mutation 



27|. 



For the same model, we investigated the state of the branching process after tobservation = 
10^ generations (i.e., time steps) for various L > 1 and a < 2 (for a > 2, the network 
turns out to be a chain). In other words, we analysed if the extinction time for given 
L and a is smaller than 10^ generations or not. On the (a/2,L/2) diagram. Fig. [TJ the 
boundary separating the extinction and non-extinction regions is a monotonously growing 
curve L{a), where L{a) oc a at sufficiently small a and diverges as a approaches 2. Note 
that L(a) actually depends on tobservation, and with increasing observation time, the area 
of extinction should increase. Based on simulations, it is practically impossible to check 
whether L(a, tobservation— ^C)o) approachs some non-trivial limit or, for any finite L and non- 
zero a, the network should finally extinct. The latter option seems to be more plausible. 

Figure [8](a) shows the evolution of Nt for a few different values of a. The averaged Nt 
(averaged over times before extinction) Nmax, decreases with a as Fig. [H](b) demonstrates. 
The simplest estimation give NmaxiO', L) ~ L/a. Figure [Hl^b) confirms that this estimate 
is reasonable, Nmax indeed inversely proportional to a, although these simulations indicate 
deviation from proportionality on L for sufficiently large L. Since new nodes are born 
uniformly randomly in the interval ( — 1, 1) from their parents, the case of L = 1 is special. 
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FIG. 6. Evolution of the number of new nodes (number of nodes in the generation t) for a single 
realization of the network of Sec. [11] defined on a one-dimensional interval —L < x < L, where 
L = l,a = 0.1. 




FIG. 7. Extinction of the network embedded in the {—L,L) ring during 10^ generations. The 
extinction and non-extinction regions are present on the L/2 vs. a/2 diagram. 
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FIG. 8. (Color online) (a) Variation of the number of nodes Nt in the current generation with 
time, for different values of a. The network is embedded in the interval —L < x < L, where L = 1, 
and only the previous generation influence the branching process. The average value of Nt at large 
t, Nmax, is represented by a solid straight line, (b) Nmax versus a for different L. 

In this situation, new nodes are actually born at any point of the ring with equal probability 
independently of the positions of their parents, and so a network structure here is not 
essential. One can consider this specific model with new nodes born in arbitrary points 
with equal probability at arbitrary L and find Nmax{0'-,L) ~ 0.5L/a. Figure [8]^b) for our 
original model shows a functionally faster growth of Nmax{0', L) with L than this proportional 
dependence. Note finally that the deviations of fluctuating A^^^ from the mean values Nmax 
in Fig. [Hl^a) are of the order of J Nmax for each L and a. 



V. NODE SPATIAL DISTRIBUTION 

In general, the nodes of the growing trees under consideration are non-uniformly dis- 
tributed in the embedding spaces. Only if the embedding area is restricted, the spatial 
distribution finally becomes uniform, see Fig. ^a). For infinite embedding space, the evolu- 
tion of the node spatial distributions is shown in Figs. [Hl^b) and (c) for the trees in which the 
birth of new nodes is determined only by a previous generation and by all existing nodes, 
respectively. The distributions in three instances are shown. The triangular shape of these 
distributions in Fig. lU^b) indicate that the spatial distribution of nodes of generation t has a 
symmetric step-function form with boarders moving away from the center (root) with con- 
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FIG. 9. (Color online) Distribution of nodes of the growing trees in space, (a) The node spatial 
distribution of the tree embedded in the interval —1 < x < +1 after 1000 time steps. The birth 
of new nodes in the tree is influenced only by a previous generation. The vertical columns of the 
histogram show the numbers of nodes within bins of width 0.2. (b) The node spatial distributions 
for the tree embedded in a one-dimensional space at different instants of the growth. The birth 
of new nodes in the tree is influenced only by a previous generation. The vertical columns of the 
histogram show the numbers of nodes within bins of width 1. (c) The same as for (b), but the 
birth of new nodes in the tree is influenced by all existing nodes. For all three plots, a = 0.01. 
Each of the results was obtained from a single realization. 



stant velocity equal approximately to 0.5, so that their coordinates increase proportionally 
to t. The density of nodes between the borders is a constant equal approximately to 0.2/a. 
In the second case. Fig. EJ^c), this expanding step-function form describes the evolution of 
the spatial distribution of all nodes in the tree. The border speed is approximately 0.6, and 
the density of nodes between borders is a constant equal approximately to 0.45/a. (Note 
that, as it should be, this value is close to the number N^ax ~ OA/a of new nodes found for 
the this tree in Sec. Illlt see Fig. |H) These observations explain the high quality of simple 
estimates obtained in Sec. [Tll 

Finally, for the networks embedded in a restricted area, in which the birth of new nodes 
is determined by a previous generation, we also measured the distribution of the number of 
nodes in one generation. We observed that this distribution is centered at Nmax and is close 
to the normal distribution. 



11 



VI. CONCLUSIONS 



We have studied the models of evolving trees embedded in a Euchdean space, in which 
the branching process is determined by the relative position of nodes in space. In these 
models, overcrowding suppresses the "fertility" of nodes. We have investigated two regimes 
of the evolution of these trees and crossover between them. In the initial stage of evolution, 
the network growth is exponentially fast, and the network is a small world. After some 
crossover time, this network becomes to grow much slowly, and, in this regime, the network 
has a large world architecture in terms of network science. We have demonstrated that 
the embedding of the network in a restricted area, which is natural for general evolution, 
sets limits to growth and can result in complete extinction. The simplest models which we 
analysed can only schematically describe real evolution processes in biology. Even these null 
models however are sufficient to demonstrate the transition from an explosive to gradual 
evolution accompanied by a dramatic change of the network structure. We believe that the 
significance of the network representation of evolutionary processes, e.g., the so-called "tree 
of life" , is greater than simply being a convenient visualization. We suggest that thorough 
exploration of the structural organization of the empirical trees of life and their analogies 
on different stages of evolution will essentially improve our understanding of evolutionary 
processes. 
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